Dissociation of quarkonium in a hot QCD medium: Modification of the 

inter-quark potential 
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We have studied the dissociation of heavy quarkonium states in a hot QCD medium by inves- 
tigating the medium modifications to a heavy quark potential. Our model shows that in-medium 
modification causes the screening of the charge in contrast to the screening of the range of the po- 
tential. We have then employed the medium-modified potential to estimate the dissociation pattern 
of the charmonium and bottomonium states and also explore how the pattern changes as we go from 
the perturbative to nonperturbative domain in the Debye mass. The results are in good agreement 
with the other current theoretical works both from the spectral function analysis and the potential 
model study. 
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INTRODUCTION 



The study of the fundamental forces among 
quarks and gluons is an essential key to the un- 
derstanding of QCD and the occurrence of different 
phases which are expected to show up when going 
from low to high temperatures and/or baryon num- 
ber densities. For instance, at small or vanishing 
temperatures quarks and gluons get confined inside 
a hadron by the strong force while at high tempera- 
tures a quite different medium consisting of quarks 
and gluons known as quark-gluon plasma (QGP) is 
expected. One of the most important features of 
the QGP formation is the color screening of static 
chromo-electric fields [l| . The suppression of heavy 
quarkonia (J/tp, Xc, ip , T) due to the color screen- 
ing analogous to Debye screening in QED plasma, 
has long been proposed as a probe of deconfinc- 
ment in a dense partonic medium. In the decon- 
fined state, the interaction between heavy quarks 
and antiquarks gets reduced due to color screening 
leading to a suppression in J/ifj yields 

Thus quarkonia at finite temperature are an im- 
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portant tool to know the status of the matter (con- 
fined / deconfined) formed in heavy ion collisions 
(see, e.g., Ref. Q). Many efforts have been de- 
voted to determine the dissociation temperatures 
of QQ states in the deconfined medium, using 
either lattice calculations of quarkonium spectral 
functions 0, [3, 0, Q or nonrelativistic calcula- 
tions based upon some effective (screened) poten- 
tials [3, [13, Ell 03 EE 0. However > the proper- 
ties of the heavy quarkonia states determined from 
the screened potentials do have a poor matching 
with the results obtained from the lattice spectral 
functions. None of the potential model studies and 
spectral functions in lattice to study quarkonia give 
a complete framework to study the properties of 
quarkonia at finite temperature. It is not yet clear 
at the moment up to what extent one may un- 
derstand the modifications of quarkonium spectral 
functions in terms of the Dcbye screening picture. 
One should not expect a precise quantitative agree- 
ment with the lattice correlators because of the un- 
certainties coming from a variety of sources [3]. 
Not only is the determination of the effective poten- 
tial still an open question but also there are other 
related issues such as relativistic effects, thermal 
width of the states and contribution from quantum 
corrections that need to be taken care of. On the 
other hand, lattice correlators are also affected by 
their own uncertainties. These may be due to the 
use of different lattices (isotropic or anisotropic). 
Additionally, the finite lattice-spacing might signif- 
icantly alter the continuum part of the spectrum. 
However, some degree of qualitative agreement had 
been found for the 5-wave correlators. This finding 
was somehow ambiguous for the P-wave correlators 
and the temperature dependence of the potential 
model was even qualitatively different from the lat- 
tice one. 
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In a recent work, Umeda [l5( found that lattice 
calculations of meson correlators at finite temper- 
ature contain a constant contribution due to the 
presence of zero modes in the spectral functions. 
The presence of a zero mode in the vector chan- 
nel had already been discussed in the literature 
while in the P-wave channels it had generally been 
overlooked. Recently Alberico et al [l6| updated 
their previous calculation [lij of quarkonium Eu- 
clidean correlators at finite temperatures in a po- 
tential model by including the effect of zero modes 
in the lattice spectral functions. These contribu- 
tions cure most of the previously observed discrep- 
ancies with lattice calculations. This observation 
supports the use of potential models at finite tem- 
perature as an important tool to complement lattice 
studies. 

The short and intermediate distance properties 
of the heavy quark interaction is important for the 
understanding of in-medium modifications of the 
heavy quark bound states. On the other hand, 
the large distance behavior of the heavy quark in- 
teraction plays a crucial role in understanding the 
bulk properties of the QCD plasma phase, viz. the 
screening property of the quark gluon plasma, the 
equation of stat e fi"7. lill and the order parameter 
(Polyakov loop) [lJJiofT 

In all of these studies, deviations from perturba- 
tive calculations and the ideal gas behavior are ex- 
pected and are indeed found at temperatures which 
are only moderately larger than the deconfinement 
temperature. This calls for quantitative nonper- 
turbative calculations. The phase transition in full 
QCD appears as a crossover rather than a 'true' 
phase transition with related singularities in ther- 
modynamic observables (in the high-temperature 
and low density regime) [21] . Therefore, it is rea- 
sonable to assume that the string tension does not 
vanish abruptly above T c . So one should study 
its effects on the behaviour of quarkonia in a hot 
QCD medium. This issue, usually overlooked in 
the literature, is certainly worth investigation. In 
the present paper, we considered this potentially 
interesting issue by correcting the full Cornell po- 
tential with a dielectric function embodying the ef- 
fects of the deconfined medium and not only its 
Coulomb part as usually done in the literature. We 
have found that this leads to a long-range Coulomb 
potential with a reduced effective charge (inversely 
proportional to the square of the Debye mass) of 
the heavy quark in addition to the usual Debye- 
screened form employed in most of the literature. 
With such an effective potential, we investigate the 
effects of different possible choices of the Debye 



mass on the dissociation temperatures of different 
quarkonium states. Since a Coulomb interaction 
always admits bound states, a criterion has to be 
adopted to define a dissociation temperature: a 
state is then considered to be melted when its bind- 
ing energy becomes of the same order as the tem- 
perature. For this purpose, we consider a gauge- 
invariant, nonperturbative form of the Debye mass 
by Kajantie et al. [22j and study systematically the 
effects of perturbative and nonperturbative terms in 
Debye mass on the dissociation pattern of quarko- 
nia in pure gauge, two-flavor, and three-flavor QCD 
respectively. Additionally, we consider the lattice 
parametrized form of the Debye mass 23 1. 



The paper is organized as follows. In-medium 
modifications to heavy quark potential is discussed 
in Section II. In subsections IIA and IIB, we study 
the medium dependence of quarkonia binding en- 
ergy and then determine their dissociation temper- 
atures in a hot QCD medium. Finally, we conclude 
in Section III. 



II. IN-MEDIUM MODIFICATIONS TO 
HEAVY-QUARK POTENTIAL 



Because of the large quark mass m = m c ,& 3> 
AqcD; the velocity of heavy quarks in the bound 
state is small and the binding effects in quarkonia 
at zero temperature can be understood in terms of 
nonrelativistic potential models [24| . More recently, 
the potential has been derived from QCD using a 
sequence of effective field theories (for a review see 
[25(1 V The present analysis also employs this idea 
to study the quarkonia with the nonrelativistic po- 
tential model. 

Let us now turn our attention to study the 
medium modifications to a heavy quark potential 
which is considered as the Cornell potential 



a 

V(r) = — 
r 



(1) 



where a and a are the phenomenological parame- 
ters. The former accounts for the effective coupling 
between a heavy quark and its antiquark and the 
latter gives the string coupling. 

The medium modification enters in the Fourier 
transform of the heavy quark potential as 



V(k) 



V(k) 
e(k) 



(2) 



where e(fc) is the dielectric permittivity given in 
terms of the static limit of the longitudinal part 
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of gluon self-energy [261 12' 

n L (o,fc,T) 



e(k) 



1 



k 2 



= 1 



" l D 
k 2 



(3) 



Note that the result for the static limit of the di- 
electric permittivity is the perturbative one. If one 
assumes nonperturbative effects such as the string 
tension survive even above the deconfinement point 
then the dependence of the dielectric function on 
the Debye mass may get modified. So there is a 
caveat about the validity of the linear dependence 
of the dielectric function (e) on the square of the 
Debye mass Mjj. For the sake of simplicity we put 
in all the nonperturbative effects together in the 
effective charge {2<j/m 2 D ) of the medium modified 
potential. 

The quantity V(k) in ([2|) is the Fourier trans- 
form (FT) of the Cornell potential. The evaluation 
of the FT of the Cornell potential is not so straight- 
forward and can be done by assuming r- as distri- 
bution (r — ► rexp(— jr)). After the evaluation of 
FT we let 7 tend to zero. Note that the FT of 
Coulomb part is straightforward to compute. The 
Fourier transform of the linear part or exp (—77*) is 



k^/2~n (7 - ifc) 3 (7 + ik) 3 



}■ 



If we put 7 = 0, we obtain the Fourier transform of 
or denoted as, 



(or) 



4a 



(4) 



Now the FT of the full Cornell potential can be 
written as 



V(fc) 



4rr 



V27rfc 4 ' 



(5) 



Substituting Eqs.(|3J) and §5§ into @ and then eval- 
uating its inverse Fourier transform one obtains the 
r-dependence of the medium modified potential [28| 
as: 



V(r) 



2a 



a) 



exp (— rnor) 



2a 



mo 



(6) 



This potential has a long range Coulombic tail in 
addition to the standard Yukawa term. In the limit 
r >> 1/tob, we can neglect the Yukawa term and 
for large values of temperature the product amp 



will be much greater than 2a/niD- So, finally the 
potential ([6]) becomes: 



V(r) 



2a 



(7) 



The above form (apart from a constant term) is a 
Coulombic type as encountered in hydrogen atom 
problem with identifying the fine structure con- 
stant e 2 with the effective charge 2a/m 2 D . Since 
mi) is an increasing function of temperature, the 
effective charge 2a/m 2 D gets waned as the temper- 
ature is increased and finally results in screening 
of the charge. The constant terms in the full po- 
tential © are introduced by hand in order to re- 
move short-distance medium effects. However, such 
terms could arise naturally from the basic compu- 
tations of real time static potential in hot QCD[29| 
and from the real and imaginary time correlators 
in a thermal QCD medium [301]. These terms in the 
potential are needed in computing the masses of the 
quarkonium states and to compare the results with 
the lattice studies. It is equally important while 
comparing our effective potential with the free en- 
ergy in lattice studies. However, these terms are 
not needed to compare the values of the dissociation 
temperatures obtained in our calculation with the 
values in lattice spectral studies because we have 
used different criteria to evaluate the dissociation 
temperatures. 

It may not be out of context to mention that the 
expression for the potential in a hot QCD medium 
is not the same as the lattice parametrized heavy 
quark free-energy in the deconfined phase (which 
is basically a screened Coulomb, for the exact form 
we refer the reader to Refs.f3ll.l32l]). As emphasized 
by Dixit [33| that one-dimensional Fourier transform 
of the Cornell potential in the medium yields the 
similar form as used in the lattice QCD to study 
the quarkonium properties which assumes the one- 
dimensional color flux tube structure. However, 
at finite temperature that may not be the case 
since the flux tube structure may expand in more 
dimensions [3lj . Therefore, it is better to consider 
the three-dimensional form of the medium modi- 
fied Cornell potential which has been done exactly 
in the present work. The medium modified poten- 
tial thus obtained has a Coulomb tail in addition 
to the screened Coulomb part. The strength of the 
Coulombic part decreases with the increase in tem- 
perature and at a certain temperature one may ig- 
nore it. 

To compare our in-medium effective potential 
with the color-singlet free-energy [3J] extracted from 
the lattice data which do not display at all any long- 
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FIG. 1: The behavior of V(r, T) /T as a function of rT for a fixed T/T c = 3.32 HI . The different curves denote the 
choice of the Debye masses. 




range term, we have plotted the full effective po- 
tential from © as a function of rT in Fig. 1. The 
lattice free energy goes to zero much faster than 
our effective potential due to the presence of the 
Coulomb tail. However, our potential employing 
the nonperturbative form of the Debye mass devi- 
ates largely from the lattice results [34| . 

Let us now proceed to study the charmonium and 



bottomonium spectrum and their binding energy 
with three possible choices of the Debye mass. Ad- 
ditionally, we take advantage of all the available 
lattice data, obtained not only in quenched QCD 
(Nf = 0) but also including two, and more recently, 
three light flavors. This enables us to study the fla- 
vor dependence of the dissociation process, a per- 
spective not yet achieved by the parallel studies of 
the spectral functions. 
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A. Binding energy of heavy quarkonia 



QCD [23J. Their forms are given below: 



Spectral function method defines binding energy 
of a quarkonium state as the distance between the 
peak position and the continuum threshold, Eu n = 
2mc,b + Voo(T) — M with M being the resonance 
mass. In our case, it is defined as the 'ionization 
potential' because of the similarity of our approx- 
imated effective potential ([7]) with the hydrogen 
atom problem. Schrodingcr equation gives the en- 
ergy eigenvalues for the ground states and the first 
excited states for charmonium (J ftp, ip' etc.) and 
bottomonium (T, T' etc.) spectra. Invoking the 
translational invariance, we can ignore the constant 
term in and the energy of the nth eigenstate is 
given by Bohr's theory: 



En 



El 

n 2 



Ei = 



rriQC 



(8) 



where vtlq is the mass of the heavy quark and Ei 
is the energy of the QQ state in the first Bohr 
state. The allowed energies for QQ states are E n = 
—Ei,—^-,—^,- ■ ■. These energies are known 
as the ionization potentials/binding energies for 
the nth bound states. It becomes a temperature- 
dependent quantity through the temperature de- 
pendence in Debye mass and it decreases with the 
increase in temperature. 

There are other states in the charmonium and 
bottomonium spectroscopy viz. x c 's and x&' s f° r 
which the determination of the medium-dependent 
binding energy is beyond the scope of our present 
calculation. For x c 's and Xb's, one should take into 
account the spin dependence of the quark-antiquark 
potential [HJ. Figures 2 and 3 show the variation 
of binding energy (in GeV) with the temperature 
(in units of critical temperatures) for J/t/j and T, 
respectively. Similar variations for other quarkonia 
(?/>' , T') can also be shown. 

In the present analysis, we consider three possible 
forms of the Debye masses, viz. the leadingorder 
term in QCD coupling {mjy), nonperturbative cor- 
rections to it (m^, p ), and the lattice parametrized 
form (m 1 ^) to study the dissociation phenomena 
of quarkonium in a hot QCD medium. The Debye 
mass at high temperature in the leading-order is 
known from long time and is perturbative [3l| in 
nature. Recently Kajantie et al. [22| computed 
the nonperturbative contributions of 0(g 2 T) and 
0(g 3 T) from a three-dimensional effective field 
theory which we consider in the present work. 
However we also consider the Debye mass obtained 
by fitting the (color-singlet) free energy in lattice 



LO 



,NP 



N . N f 



Ng 2 T 

47T 



In- 



g 2 T 



' D 



= I Am 



to 
D 



c N g 2 T + d NtNf g 3 T 
(9) 



where the coefficient c N captures the nonperturba- 
tive effects and d,N,N f is related to the choice of the 
scale in rn^°. We employ the two-loop expression 
for QCD coupling constant at finite temperature [36j 
and choose the renormalization scale determined in 

Note that different curves in each figure denote 
the choice of the Debye masses (|9]) used to calcu- 
late the binding energy from ((8]). There is a com- 
mon observation in Figs. 2,3 that there is a strong 
decrease in binding energy with the increase in tem- 
perature. In particular, binding energies obtained 



from 



and nip give realistic variation with the 



temperature. The temperature dependence of the 
binding energies shows a qualitative agreement with 
the similar variations shown in [23| . 

However, when we employ the nonperturbative 
form of the Debye mass rn^ p the binding energies 
become unrcalistically small compared to the bind- 
ing energy at T = and also compared to the bind- 
ing energies employing ni£° and m^,. This can be 
understood by the fact that the value of m^ p is 
significantly larger than both mfj and m^. This 
observation hints that the present form of the non- 
perturbative corrections to the Debye mass may not 
be the complete one, the situation may change when 
nonperturbative corrections of higher order 0(g 4 T) 
are added to the Debye mass and then use it to 
calculate the binding energy. 

Thus a study of the temperature dependence of 
the binding energy is poised to provide a wealth 
of information about the dissociation pattern of 
quarkonium states in thermal medium which will 
now be used to determine the dissociation temper- 
atures of different states. 



B. The dissociation temperatures for heavy 
quarkonia 



Dissociation of a two-body bound state in a ther- 
mal medium can be understood as: when the bind- 
ing energy of a resonance state drops below the 




mean thermal energy of a parton, the state be- 
comes feebly bound. The thermal fluctuations then 
can destroy it by transferring energy and exciting 
the quark-antiquark pair into its continuum. The 
spectral function technique in potential models de- 
fines the dissociation temperature as the tempera- 
ture above which the quarkonium spectral function 
shows no rcsonancc-likc structures but the widths 
shown in spectral functions from current potential 
model calculations are not physical. The broaden- 
ing of states with the increase in temperature is not 
included in any of these models. In Ref.[23[, the au- 
thors argued that one need not to reach binding en- 
ergy (Ebin) to be zero for the dissociation. Rather 
a weaker condition Eun < T causes a state weakly 
bound and the thermal fluctuations can then de- 
stroy it. Since the (relativistic) thermal energy of 
the partons is 3T, the lower bound on the dissocia- 
tion temperature (Tjj) is obtained from the relation 

V^=3Ib, (10) 

where string tension (a) is taken as 0.184 GeV 2 , 
and critical temperatures (T c ) are taken 27QM eV, 
203MeV and 197 MeV for pure, two-flavor and 
three-flavor QCD medium, respectively [1^. The 
dissociation temperatures for the ground states and 
the first excited states of cc and bb are listed in 
Table I with the Debye mass in the leading-order 
m L °. It is seen from Table I that ip' is dissoci- 
ated in the vicinity of critical temperature while 
J/tp and T' are dissociated around 1.2 T c . T is dis- 
sociated at a relatively higher temperature 1.6T C . 
These values agree quantitatively with the recent 
values reported by Mocsy and Petreczky (23| . On 
the other hand, when we use Debye mass from lat- 



TABLE I: Lower (upper) bound on the dissociation 
temperature(r_D) for the quarkonia states (in units of 
T c ) using the leading-order term in the Debye mass 



State 


Pure QCD 


N f = 2 


N f = 3 


JH> 


1.1 (1.4) 


1.3 (1.7) 


1.2 (1.6) 


</>' 


0.8 (1.0) 


0.9 (1.2) 


0.9 (1.1) 


T 


1.4 (1.8) 


1.7 (2.1) 


1.6 (2.0) 




1.0 (1.3) 


1.2 (1.6) 


1.2 (1.5) 



TABLE II: same as Table I but using full potential (0 . 



State 


Pure QCD 


N,=2 


N f = 3 


J/V 


1.2 (1.5) 


1.3 (1.72) 


1.3 (1.7) 




0.8 (1.2) 


1.0 (1.4) 


1.1 (1.2) 


T 


1.4 (1.8) 


1.7 (2.3) 


1.6 (2.1) 




1.0 (1.4) 


1.2 (1.7) 


1.2 (1.5) 



tice parametrized free energy (nip), the values be- 
come much lowered than the leading-order results 
(see Table III). However, nonperturbative correc- 
tions to the Debye mass (m^ ) p ) make the values 
unrealistically small. These observations can be 
understood from the hierarchy in their numerical 
values: < < m^ p . 

The fact that m^ p leads to unrealistic smaller 
values of dissociation temperatures does not imply 
that one should ignore the nonperturbative terms 
in the Debye mass. In fact, nonperturbative terms 
cannot be ignored in the regime where coupling is 
strong which is indeed the case dealt with. It would 
rather be of interest to raise the question why this 
nonperturbative result obtained with a dimensional 
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TABLE III: Lower (upper) bound on the dissociation 
temperatures using the lattice parametrized form of the 
Debye mass m D . 



State 


Pure QCD 


N f =2 


N f =3 


J/i> 


0.8 (1.0) 


0.9 (1.2) 


0.9 (1.1) 


$ 


0.5 (0.7) 


0.7 (0.8) 


0.6 (0.8) 


T 


1.0 (1.3) 


1.2 (1.6) 


1.2 (1.5) 




0.7 (0.9) 


0.9 (1.1) 


0.8 (1.0) 



reduction is not in agreement with the Debye mass 
arising from Polyakov-loop correlators. This could 
be partially due to the arbitrariness in the definition 
of dissociation temperature, since strictly speaking 
a Coulomb potential always admits bound states 
in its spectrum. Indeed the choice of the average 
thermal energy 3T is not rigid because even at low 
temperatures T < T c (say) the Bose/Fermi distri- 
butions of partons will have a high energy tail with 
partons of mechanical energy greater than the bind- 
ing energy. So, we calculate the upper bound of the 
dissociation temperatures by replacing the average 
thermal energy ~ T which is listed within the first 
bracket in the tables where the values are increased 
by 30% approximately. 

The results for the dissociation temperatures of 
various quarkonia listed above in Table I and Ta- 
ble III are obtained by dropping all the finite-range 
terms in the full effective potential ([6]). As men- 
tioned earlier, for the s-wave states this leads to an 
analytically solvable Coulomb potential. To see the 
effects of the finite-range terms in ©, we solve the 
Schrodinger equation numerically with the full ef- 
fective potential© and determine the energy spec- 
trum of the ground and the first excited states of 
the charmonium and bottomonium spectrum with 
the Debye mass in the leading-order. We find that 
the dissociation temperatures change by ~ 10 20% 
shown in Table II. The dissociation temperature for 
J/^> in the pure gauge case becomes 1.2T C which 
was earlier 1.1 T c (Table I) and for tp', it now be- 
comes 1.0T C for Nf = 2 which was earlier 0.9T C . 
The same trend follows for other charmonium and 
bottomonium states. This slight increase in the dis- 
sociation temperatures is caused by the increase in 
the binding energies due to the finite-range terms 
in the full potential 

Finally, to compare our results with a recent cal- 
culation [23j having the same input based on a po- 
tential study for a three-flavor QCD with T c = 
192MeV, we calculated the upper bound of disso- 
ciation temperatures in Table IV and V with the 
same form of Debye mass used in Ref. [23| . It shows 



TABLE IV: Upper bound on the dissociation temper- 
atures (T D ) with T c = 192 MeV [U using the lattice 
parametrized form of the screening mass(m^). 



State 


i>' 


J/V> 




T 


To 


< 0.9T C 


1.2T C 


1.1T C 


1.6T C 



TABLE V: Same as Table IV but using the full potential 



State 


i>' 


J/i> 




T 


T D 


< 1.0T C 


1.3T C 


1.2T C 


1.6T C 



a good agreement with their results [23j. However, 
the agreement holds well even with the full Cornell 
potential. 



III. CONCLUSIONS AND OUTLOOK 

In conclusion, we have studied the dissociation 
phenomena of quarkonia in a hot QCD medium 
by investigating the in-medium modifications to a 
heavy quark potential. We have found that medium 
modification causes a dynamical screening of color 
charge which in turn, leads to a temperature de- 
pendent binding energy. We have systematically 
studied the temperature dependence of binding en- 
ergy for the ground and first excited states of char- 
monium and bottomonium spectra in pure and re- 
alistic QCD medium. We have then determined 
the dissociation temperatures employing the per- 
turbative result of Debye mass (m)^) and the lat- 
tice parametrized form m 1 ^. Our estimates are con- 
sistent with the finding of recent theoretical works 
based on potential models [23| . However, these val- 
ues are signifi cantly smaller than the predictions of 
Refs. pi [tiH. [iol] based on the first principle lattice 
calculations which are however plagued by its inher- 
ent uncertainties. In contrast, the inclusion of non- 
perturbative contributions to Debye mass lowers 
the dissociation temperatures substantially which 
looks unfeasible. Thus, this study provides us a 
handle to decipher the extent up to which and how 
much nonperturbative effects should be incorpo- 
rated into the Debye mass. 

In brief, J ftp is found to be dissociated at temper- 
ature above the critical temperature (around 1.2T C ) 
when the lcadingordcr term in the Debye mass has 
been employed. However, it is dissociated just be- 
low the T c when lattice parametrized form of Debye 
mass (nonperturbative) has been employed. This 
finding ensues a basic question about the nature of 
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dissociation of quarkonium in a hot QCD medium. 

Finally, our approach based on the in-medium 
modifications provides charmonium and bottomo- 
nium dissociation temperatures which agree nicely 
with recent quarkonium spectral function studies 
using a potential model |23[. This is true only for 
the perturbative result for the Debye mass but non- 
perturbative corrections to it make the melting tem- 
peratures too low to compare to the spectral analy- 
sis of the lattice temporal correlator of the mesonic 
current. This leaves an open problem of the agree- 
ment between these two kinds of approaches. This 
could partially be due to the arbitrariness in the 
definition of dissociation temperature. To examine 
this point we have estimated the upper bound of 
the dissociation temperatures from the condition: 
Ebin = T. We found that these estimates obtained 



by employing the lattice parametrized Debye mass 
show good agreement with the predictions in (23j 
which was not true for the earlier definition (| 1 01) : 
Ebin = 3T. However, a numerical solution of the 
Schrodingcr equation with the full effective poten- 
tial ([6]) gives in general a slightly higher value of the 
dissociation temperatures. For Xc(Xfc) melting tem- 
peratures, one would start with the spin-dependent 
heavy quark potential and follow the same proce- 
dure. We will look into it in the future. Finally 
it would be of interest to study the corresponding 
quarkonium spectral function (and temporal corre- 
lator) after giving these states a thermal width. 
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